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Abstract. We study the free energy of a heavy quark-antiquark pair in a thermal medium. We constuct a 
simple ansatz for the free energy for two quark flavors motivated by the Debye-Hiickel theory of screening. 
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1 Introduction 

The theory of quantum chromodynamics (QCD) predicts 
that at high temperature and density there is phase tran- 
sition from hadron gas to quark gluon plasma (QGP). 
This QGP phase is believed to be formed in the colli- 
sion of very high energy heavy-ions. Various probes have 
been proposed to detect the formation of QGP in high 
energy heavy-ion collision experiments. One of the clean 
signatures of QGP formation is the suppression of the 
J/W P. At high temperatures and density, medium ef- 
fects screen the interaction between the heavy quark and 
antiquark pairs (QQ). The Debye screening radius charac- 
terizes the distance beyond which significant thermal mod- 
ification arise in the heavy quark potential. Qualitatively, 
one may expect that the QQ dissociates in the medium 
when the Debye radius becomes smaller than the radius of 
QQ bound state. Quantitative predictions, such as disso- 
ciation temperatures or energy densities, however, require 
the study of quarkonia in a medium. Recently there have 
been some finite temperature studies of QQ spectra in 
quenched lattice QCD [2] . Previously, the medium depen- 
dence of quarkonia had been studied in the Schrodinger 
equation formalism. For these calculations, one needs the 
potential between the heavy quark and antiquark in the 
medium. However, this potential is not yet well under- 
stood. Only for very small distances or very high temper- 
atures one can calculate it perturbatively. On the other 
hand, the free energy F(r,T) of a static quark-antiquark 
pair can be calculated in lattice QCD, and using thermo- 
dynamic relations, one can then derive the change in en- 
ergy of a thermal medium due to the presence of a QQ pair 
from this free energy. This may be more closely related to 
the desired potential. In any case, it is desirable to under- 
stand how the free energy F(r, T) depends on r and T. So 
far, a functional form of F(r,T) in accord with the lat- 
tice data over wide range of temperatures is still missing. 
In this work we construct such a form of F(r,T), based 



on studies of screening in the Debye-Hiickel theory ,5 . Al- 
though there are only two parameters in our analytical 
form of F(r, T), we can fit the lattice data quite well for 
all r and for a wide range of temperatures, < T < 2 T c . 

This paper is organized as follows. In section 2 we 
discuss the screening of the free energies of static quark 
anti-quark pair in the presence of medium and derive our 
ansatz for F(r,T). In section 3 we fit our form of free 
energy F(r, T) to the lattice results, thus obtaining the 
temperature dependence of the free parameters in F(r, T). 
Section 4 contains discussions and conclusions. 



2 Medium effects on the interaction of QQ 

The free energy of a static quark-antiquark at zero tem- 
perature is given by [3] 



F(r,T = 0) = err - -, 
r 



(1) 



where a is the string tension and a is the gauge coupling 
constant. The first term is due to formation of flux tube 
or string between the Q and the Q when they are pulled 
apart. This string free energy increases linearly with the 
separation of the QQ. The second term in the above equa- 
tion corresponds to the Coulomb interaction between the 
static charges. This free energy F(r,T = 0) describes 
the charmonium and bottomonium spectrum quite well 
in terms of a and a. 

The free energy F(r, T) increases linearly with r un- 
til at some ro, it becomes energetically favorable for the 
string to break into two heavy-light mesons (D, B) rather 
than stretching further. After the string breaks at ro, the 
free energy remains constant for all separations r larger 
than ro. The value the free energy F(ro,T — 0) at the 
string breaking point ro can be calculated. Thus, the masses 
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of the J/\P, the D and F(r Q ,T = 0) are related by 



2M 



D 



M 



7/f 



EB 



= 2m c + F(r Q ,T = 0). 



(2) 



Here EB denotes the binding energy and Md, Mj/y and 
m c the masses of the D, the J/<F and the charm quark, 
respectively. One can write a similar equation for the bot- 
tomonium system. Rearranging the terms, one arrives at 
the following relation 



F(r ,T = 0) ~ 2{M D - m c ) ~ 2{M B - m b ). 



(3) 



Assuming universal behavior of the string one will get a 
universal value of rg, as both 2(Mrj — m c ) and 2(Afg — mj) 
are roughly equal to 1.2 GeV. Using the conventional value 
of cr ~ (0.4 GeV) 2 for the string tension, we get 



ro 



2(M D , B - m C)6 ) 1.2GeV 



1.5fm 



At finite temperature, medium effects screen the heavy 
quark-antiquark interactions. The mechanism is quite well 
understood for Coulombic systems and leads to the well- 
known Debye-screened Coulomb potential and the corre- 
sponding Coulombic free energies @]. In general we may 
separate the free energy F(r, T) of a pair of static charges 
into contributions arising from its potential energy at zero 
temperature F(r, T = 0), and a screening function f(r, T), 



F(r,T)=F(r,0)f(r,T) 



(4) 



This can be formulated in the framework of Debye-Huckel 
theory 0] and can easily be solved for Coulombic sys- 
tems. The Debye-Huckel ansatz has been generalized to a 
large class of potentials in arbitrary dimensions [5] . In the 
following we briefly discuss the formalism of the Debye- 
Huckel theory and will then come back to the problem 
of finding an appropriate screening ansatz for the heavy 
quark potential. 



2.1 Debye-Huckel Theory 

In the framework of Debye-Hiickel theory, the free energy 
F(r, T) can be derived from a corresponding Poisson equa- 
tion 0. For example, a free energy F at T = with the 
form, 



F = (3^ 



(5) 



in three space dimensions can be obtained by solving the 
Poisson equation 



V 2 F (v + l) 



rV +2 



VF ■ f = 47r/3<5(r). 



(6) 



In the medium, the interaction between the heavy quark- 
antiquark changes because the heavy quarks polarize the 



medium. In linear response theory this leads to the follow- 
ing substitution for the source/charge term in the Poisson 
equation 



5(r) — ► 5(r) + AF. 



(7) 



All effects of the medium are contained in the factor A. In 
the case of abelian charges, A is function of the density of 
charges in the medium. In a non-abelian medium, in which 
we are interested, making an estimate of the density of 
color charges is very difficult. However, if the free energy 
is known, one can give an estimate of A. This is what we 
basically do in this work. A has the mass dimension of 2 in 
3-dimensional physical space. Defining the screening mass 
[i = (AirftA) 1 ' one arrives at the following modified 
Poisson equation, 



1 d 2 F l-rjdF 



(8) 



The solution of this equation gives the free energy of static 
charges in a thermal medium interacting through a poten- 
tial described by equation (5). The equations in terms of 
screening functions can be obtained by substituting equa- 
tion Q in the above equations. In the following we con- 
sider these separately for the case of string and Coulomb 
free energies, by choosing the appropriate values of /?, 77 in 
equation (5). 



2.2 Coulomb Screening 

In this section we obtain the r and T dependence of the 
Coulomb screening function f c . The temperature depen- 
dence appears only through the screening mass n, as one 
can see from equation (8). For Coulomb interactions, the 
two parameters of equation (5) are (3 = —a, rj — — 1. 
Substituting the form F c (r,T) = -(a/r)f c (r,T) for the 
Coulomb free energy in equation (8), we get the following 
equation 



ld 2 f c 



dr 2 



fc 



-4irS(r) 



(9) 



for / c (r, T). It is easy to see from this equation that f(r, fi(T)) 
will depend only on the combination /ir. The solution is 
given in terms of the modified Bessel function Ki/ 2 and 
becomes by, 



fc 



2 i/a 



HrKi (/if) 



(10) 



which is in the familiar Debye screening form. However, in 
the presence of medium the Coulomb free energy has to 
satisfy an additional boundary condition when the static 
charges are infinitely far apart; it must then approach the 
value F[r — 00, T) = —afi for each pair of opposite color 
charges |1]E]. Taking this fact into account, the Coulomb 
screening function is given by 



fc = e 



\ir. 



(11) 
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This form satisfies all the required boundary conditions. 
The Coulomb free energy in the presence of medium is 
thus given by 



F c (r,T) 



a r _ 

-- e 
r L 



■fxr 



(12) 



2.3 String Screening 

Next we will derive the string screening function. As in 
the Coulomb screening case, the temperature dependence 
enters in the free energy only through the screening mass 
(i. For the string interaction, the parameter (3 of equation 
(5) is the string tension a, and 77 = 1. We substitute the 
form F s (r, T) — arf s (r, T) for the string free energy in 
equation (8), which gives the following Poisson equation 
for f s (r,T), 



I d 2 /. 
r dr 2 



~?-T-~ /j 4 r/ s = 47T<5(r) 



(13) 



Here, as in the case of equation (9), the r and \x depen- 
dence are such that the solution f s is a function of the 
product /ir only. The solution of equation (13), which sat- 
isfies the boundary condition f(r, /i = 0) = 1, is given 

by 



fs{x) 



HI/4) ^ 
2 3 / 2 r(3/4) 2 3 /4r(3/4) 



#1/4 (x 2 ) 



,(14) 



with x = (ir, and where -K1/4 is the modified Bessel func- 
tion. The string free energy is now given by 



F s (r,T) = 



£0/4) ^ 
2 3 / 2 r(3/4) 2 3 / 4 r(3/4) 



K 1/A (x 2 ) 



.(15) 



Because of the medium effects, F s (r,T) eventually stops 
rising linearly with r and flattens off for large r, giving 



F s (r = 00,(1) = 



* £(1/4) 
/j2 3 /2r(3/4) 



(16) 



It is important to note that the treatment of the medium 
effects in equation (7) is an approximation. In principle, 
there can be higher order corrections to equation (7) if the 
medium is not dilute enough. Also one may expect non- 
trivial confinement effects on equation (7). As a simple 
ansatz, we here assume that all such effects only change 
the argument of the modified Bessel function K 1 u, 



#1/4 (x 2 ) -> K 1/4 (x 2 + kx 4 ) 



(17) 



The resulting higher order corrections affect the solution 
in the intermediate range of x = (ir. The string free energy 
in the medium now takes the form, 



F(r,T) 



^(1/4) 
23/2/^(3/4) 2 3 /4r(3/4) 



K 1/4 (x 



(18) 



2.4 A Model for Screening of the QQ Free Energy 

In the above two subsections we have discussed separately 
the screening of the two terms contributing to the heavy 
quark potential given in equation^ In general, we should 
now solve the Poisson equation arising in the Debye-Huckel 
theory for the combined potential. This, however, seems to 
be difficult to achieve in closed form. At present, we there- 
fore explore a screening ansatz for the heavy quark free 
energies based on the asymptotic forms discussed above, 
i.e., we assume that both terms are modified by their ap- 
propriate screening functions. 

This leads to the modification of the free energy F(r, T) 
of the form 



F(r,T) = arf s (r,T)--f c (r,T). 

r 



(19) 



This free energy is same as the zero temperature form, 
modified by the Coulomb and string screening functions, 
respectively. These screening functions satisfy the bound- 
ary conditions 



fs(r,T) 
fs(r,T) 



fc(r,T) = l, 
/ c (r,T) = l, 



r -> 0. 
T ^ 0. 



(20) 



The first condition arises because for distances r << 1/T 
the QQ cannot see the medium, the second because there 
is no medium at T = 0. 

The total free energy in the medium is then given by 
the sum of the screened Coulomb and string free energies. 
Using equation (12) and (18) we get 



F(r,T) 



r(i/4) 



2 3 /2r(3/4) 23/4^(3/4) 
[e- x + x] . 



K 1/4 (x 2 



(21) 



In the following we will compare this ansatz with lat- 
tice results for the free energies of a static quark anti-quark 
pair in a thermal heat bath of quarks and gluons. 



3 Numerical Results 

The free energy F(r, T) given in equation (|21|l contains 
only two temperature dependent parameters, (i and k; we 
assume the string tension a and the gauge coupling a to be 
temperature-independent. Since we are not analyzing the 
heavy quark free energies at very short distances, where 
the running of the coupling becomes relevant, we use for 
a the string model value a = n/12. To determine the 
temperature dependence of (i and k, we fit F(r,T) to the 
lattice results obtained in 2-flavor QCD 7 . At r = 00, the 
free energy F(r, T) depends only on fi and not on k, 



F(T) 



r(i/4) 



H(T) 23/2/^(3/4) 



a(x(T) 



(22) 



with F(r = 00, T) = F(T). To obtain the temperature 
dependence of (i, we fit this F(T) to the lattice data for 
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Fig. 1. Results from 2-flavor lattice QCD calculations for the 
large distance limit of the free energy (upper figure) and the 
corresponding screening mass (lower figure) 
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Fig. 2. Fits to the lattice results for the free energy as func- 
tion of r for different temperatures (upper figure), and for the 
parameter k (lower figure). 



r = oo. For /u(T), this gives the following form as function 
0fF(T), 



+ V / ( f1(t)+4 '5#w3)) /4q2 ' ,23) 

In Fig. 1 we show the lattice results for F(T), obtained 
from the mentioned lattice calculations in 2-flavor QCD 
\7\ , as well as the corresponding fit results for the screening 
mass jtt(T). 

Once we have the temperature dependence of fJ-(T), we 
fit equation (21) to the lattice data to obtain k(T). Our 
form of the free energy F(r, T) fits these data quite well, as 
seen in Fig. 2, where we show our fit curves (solid lines) 
together with the lattice results. The fit results for the 
temperature dependence of n(T) are also shown in Fig. 
2. They indicate that k(T) is finite and nearly constant 
below T < 0.8T C , which may reflect the effect of confine- 
ment and of string breaking as an additional non-thermal 
screening effect, which we have ignored here. Closer to T c , 
k(T) drops sharply and remains small above T c . In fact, 
above T c we obtain good fits also with k set to zero (as 
shown by the lower lines for T=T C and 1.02 T c in Fig. 3). 
This suggests that our ansatz with two separate screen- 
ing functions provides a good approximation to the correct 
screening solution of the Debye-Hiickel theory in the QCD 
plasma phase. 



4 Discussions and Conclusions 

Using Debye-Hiickel theory, we have constructed a sim- 
ple functional form for the free energy F(r, T) of a static 
quark-antiquark pair in a thermal heat bath. For con- 
stant string tension and gauge coupling constant, F(r,T) 
is specified by only two temperature-dependent parame- 
ters, n and k. We determine their temperature dependence 
by fitting our form of F(r, T) to the lattice results for the 
free energy. The resulting F(r, T) fits the lattice data quite 
well for all r and in a broad range of temperatures from 
zero to 2 T c . Since our F(r,T) is an analytic form, one 
can now obtain from it the potential energy of a static 
quark anti-quark system using the thermodynamic rela- 
tions. This may then be applied to study medium effects 
on the properties of heavy quarkonia. 
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